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POPULATION SYNTHESIS FOR SYMBIOTIC STARS 
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ABSTRACT 

We have carried out a detailed study of symbiotic stars with white dwarf accretors by 
means of a population synthesis code. We estimate the total number of symbiotic stars 
with white dwarf accretors in the Galaxy as 1,200 - 15,000. This range is compatible 
with observational estimates. Two crucial physical parameters that define the birthrate 
and number of symbiotic stars are the efficiency of accretion by white dwarfs (which 
greatly depends on the separation of components after common envelope stage and 
stellar wind velocity) and the mass of the hydrogen layer which the white dwarf can 
accumulate prior to the hydrogen ignition. The theoretical estimate of the Galactic 
occurrence rate of symbiotic novae ranges from about 1.3 to about 13.5 yr -1 , out of 
which weak symbiotic novae comprise about 0.5 to 6.0 yr -1 , depending on the model 
assumptions. We simulate the distributions of symbiotic stars over orbital periods, 
masses of components, mass-loss rates of cool components, mass-accretion rates of hot 
components and luminosity of components. Agreement with observations is reasonable. 
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1 INTRODUCTION 

Symbiotic stars (SySs) are an inhomogeneous group of vari- 
able stars with composite spectra. Their spectacular spec- 
tral and photometric variability is a very important and in- 
teresting phenomenon. The spectra of SySs suggest that a 
three-component system consists of a binary system in which 
an evolved giant transfers matter to a much hotter com- 
pact companion by means of stellar w ind and an HII region 
i|BermarJll93^ : lBovarchuklll967l. ll96Sl) . 

The cool component is a red giant (RG) which 
is a first giant branch (FGB) or an asymptotic giant 
branch (AGB) star. In the majority of SySs the hot 
component is, most probably, a white dwarf (WD), a 
sub d warf or an accreting low-mass main-sequence (MS) 
star llTutukoy fc Yungelsonlll976tlKenvon fe Webbinkll983 



Miirset ct al. 1991; Yungelson ct. al. 1995, hereafter YLTK; 



Iben fe Tutukovlll996l) . The variability of SySs may be due 



to th e thermonuclear runaways on t he surface of an accreting 
WD (iTutukov fc Yungelsonl Il976t lYLTKt llben fc Tutukol 



1996) or to the variations i n the accretion 
the hot component llDuschll Il986l 



rate onto 
Bisikalo et"afl [2002 



iMitsumoto et all 120051) . Recent reviews of the properties 
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of SySs can be found in IMiirset IT Schmidl Jl999i) and 
iMikolaiewskal J2003ft . The nova-like eruptio ns of SySs were 
reviewed bv lMikolaiewska fc Kenvonl (Il992t) . 

In astrophysics today, the interaction between compo- 
nents in the binaries is of s pecial interest , for which SySs of- 
fer an exciting laboratory llKenvon1ll986h . For instance, the 
nature of components' interaction and of the activity of WD 
accretors are still controversial (IMiirset fc Schmidl 11999): 
mass loss and stellar wind (including the terminal velocity 
and acceleration mechanism) from a cool giant are not clear; 
SySs have b een frequently discussed as Type la supernovae 
progenitors (Tutukov & Yungelson 1976; iMunari fc Renzinil 
I l992t iKenvon et al l Il993t lYLTKt llben fc Tutukovl Il996t 
lHachisu et al.l ll999). but their contribution to the rate of 
these events remains uncertain. 

Theoretical studies on the formation and e volution of 
the Sy S s have been pub l ished , e.g., by YLTK. [Han et ail 
<1995a^liben fc Tutukol Jl99rJ) : lHurlev et alJ (120021) . Their 
investigations reproduced successfully many observed prop- 
erties of these objects. However, in the recent years, a num- 
ber of papers have appeared that provide new insights in 
the different aspects of st ellar evolution relevant t o the SySs. 
iPrialnik fc Kovetzl dl995l) and lYaron et all (l2005Tt did a more 
detailed study of Novae m odelsT lNel^m^nser^n d2000f) and 
iNelemans fc Tout] ([2005) put forward an alternative algo- 
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rithm for angular momentum loss during the common- 
envelope stages of evolution; Wi nters and his collaborators 
ijWinters et al.ll2000l 120021 l2003h studied the hydrodynami- 
cal structure of the stellar wind around AGB stars with low 
mass-loss rate and low wind outflow velocity. New observa- 
tional catalogue and an ana lysis of SySs have be en given in 
iBelczvnski et alJ l l200Ch and lMikolaiewskal J2003I) . It is time 
for a new study of the SySs. 

In the present work we model the subpopulation of SySs 
with WDs as the hot components and FGB or AGB stars as 
the cool components. We study evolutionary channels lead- 
ing to the formation of symbiotic systems in which He, CO 
or ONe WDs accrete hydrogen-rich matter from the stel- 
lar wind of FGB or AGB stars. We obtain SySs' birthrate, 
lifetime and number in the Galaxy, and some potentially ob- 
servable parameters of SySs, such as orbital periods, masses 
of the components, their luminosities and mass-loss rates. 
Special attention is paid to the dependence of population 
model on the parameters entering population synthesis. This 
paper is the first one of a series of papers on SySs. It offers 
a basis for further studies of symbiotic stars, such as D-SySs 
which are SySs with thick dust shells (paper in preparation) . 

In § 2 we present our assumptions and describe some 
details of the modeling algorithm. In § 3 we discuss the main 
results and the effects of different parameters. In § 4 the 
main conclusions are given. 



2 MODELS 

In SySs, the cool component loses matter at a high rate by 
stellar wind and the hot component moves in the wind and 
accretes enough matter to produce symbiotic phenomenon. 
In this paper, binaries are considered as SySs if they sat- 
isfy the following conditions: (i) The systems are detached; 
(ii) The luminosity of the hot component is larger than 
lOL© which is the "thresho ld" luminosity for the hot compo- 
nent of SySs a s infe rred by Mikoi alewska fc Kenvonl il992l) , 
iMiirset et all Jl99lf) and YLTK; it may be due to the ther- 
monuclear burning (including novae outbursts, stationary 
burning and post-eruption burning) or the liberation of grav- 
itational energy by the accreted matter; (iii) The hot com- 
ponent is a WD and the cool component is a FGB or an 
AGB star. Below we describe our selection algorithm and 
give some details that are important for understanding the 
model. All estimates made in the paper are for stars with 
X=0.7 and Z=0.02. Chandrasekhar mass limit is taken as 
1.44 M . 
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Figure 1. — Evolutionary channels for formation of SySs. From 
left to right channels I, II and III, are shown. RLOF stands for 
the Roche lobe overflow, MS for main sequence, WD for white 
dwarf, FGB for the first giant branch, AGB for asymptotic giant 
branch, SySs for symbiotic stars. Numbers represent fractions of 
the initial sample of binaries that proceed through the consecutive 
evolutionary stages of each channel for the "standard" case 1 (see 
Table 0. 



channel, the formation of a helium star is possible. If the he- 
lium star is hot enough and the matter inflow rate into the 
circum-binary medium is high enough, then symbiotic phe- 
nomenon can arise. However, the number of systems with 
helium stars is very sm all compared to the numbe r of sys- 
tems with white dwarfs jYLTKt iHurlev et afll2002ft and we 
omit them from consideration. 

Through channel III go the systems that are initially 
wide. In these systems only CO WD or ONe WD accretors 
are formed. 

We consider ONe WDs as remnants of AGB stars that 
still experience non-degenerate carbon ignition, but avoid 
elect ron captures on N e and Mg in the core. We assume 
after iPols et all jl998T) that for solar metallicity stars the 
corresponding range of masses is 6.1 - 7.9 Mq. However, 
the range of masses of ONe WD progenitors still remains 
uncertain and it may be shifted to highe r masses and be 
much more narrow, several 0.1 Mn o nly ( Ibgji_fcjrutukov 



1985; 
2003; 



Gil-Pons fc Garci'a-Berrol I2OOIL I2OO2I iGil-Pons et 



Siessfl^od) . It is possible that we overestimate the 



number of systems with ONe WD and, more g enerally, 
initial-final mass relation of P ols et alJ (Il998h is too 
steep. 



2.1 Formation Channels of Symbiotic Stars 

All progenitors of SySs pass through one of the three routes 
(YLTK): (i) unstable Roche lobe overflow (RLOF) with for- 
mation of a common envelope; (ii) stable RLOF; (iii) for- 
mation of a white dwarf+giant pair without RLOF. These 
scenarios are shown in Fig. 

Through channel I pass the systems with short orbital 
periods. The primary overflows its Roche lobe in the FGB 
or AGB stage and forms a common envelope. After ejecting 
the common envelope, it transforms into a WD. 

In channel II, the systems undergo stable RLOF. In this 



2.2 Common Envelope Evolution 

In channels I and II, the primary can overflow its Roche lobe. 
If the mass ratio of the components (q — Maonor/Afaccretor) 
at the onset of RLOF is larger than a certain critical value 
g c , the mass transfer is dynamically unstable and results 
in the formation of a common envelope. The issue of the 
criterion for dynamically un stable RLOF q c is still open. 
iHiellming fc Webbinkl ( 1987) did a detailed st udy of sta- 
bility of ma ss transfer using polytropic models. lHan et al] 
feOOll |2002|] showed that q c depends heavily on the as- 
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sum ed mass-trans f er effi ciency. In this work, we adopt 
after iHurlev et alJ <2002l) : 



<7c 



1.67 -a + 2 



Mc 
M 



/2.13, 



(1) 



where M c and M are core mass and total mass of the donor, 
respectively and x = 0.3 is the exponent of th e mass-radius 
relati on at constant luminosity for giant stars ijHurlev et alJ 
l2000l) . 

For the common envelope evolution, it is generally as- 
sumed that the orbital energy of the binary is used to expel 
the envelope of the donor with an efficiency a cc : 



eA-Eori: 



(2) 



where i?bind is the total binding energy of the envelope and 
A_E or b is the orbital energy released in the spiral-in. In the 
present paper w e apply Eq. in the form suggested by 
IWebbinkl Jl984) with modifications after Ide KoolT(ll990l) : 



G(M C + M e )M e _ ( GM c m 



GMm 



(3) 

Here A is a structure parameter that depends on the evo- 
lutionary stage of the donor, ai is the orbital separation at 
the onset of the common envelope, M, M c , and M B are the 
masses of the donor, donor's envelope and the core, respec- 
tively, Ri is donor's radius and m is the companion mass. 
Then the orbital separation of a binary after common enve- 
lope phase at is given by 



Mc 
M 



1 + 



2Af c ai 



(4) 



a ce Am.Ri 

iNelemans et alJ ll2000F) suggested to describe the vari- 
ation of the separation of components in the common en- 
velopes by an algorithm based on the equation for the system 
orbital angular momentum balance which implicitly assumes 
the conservation of energy: 

AJ M e 



: 7- 



(•») 



M + m' 

where J is the total angular momentum and AJ is the 
change of the total angular momentum during common en- 
velope phase. The orbital separation ai after common enve- 
lope phase is then 

2 



(j. f 



M \ ( M c + 



M c I \ M 



1-7 



M c 



M + m 



(6) 



In the runs computed with Eq. © we apply it irrespec- 
tive of the mass ratio of components, as suggested by 
INelemans fc Touti " 



Following INelemans fc Toutl ( 2005) we call the formal- 
ism of Eq. © a-algorithm and that of Eq. 7-algorithm. 
For a-algorithm, there are two parameters: a ce and A. Both 
parameters are highly uncertain. It's still not completely 
clear whether sources other than gravitational energy have 
to be taken into account when computing a co and A and 
how the core-envelope boundary for the estimate of A has to 
be defined (see, e. g., lHan et al.1 1995b; Ibcn & Livio 1993; 
[ Dewi fc Tauris l200fllTauris fc Dewi l200lHSoker fc Haroal 
2003). Both parameters depend on the evolutionary stage 
in which RLOF occurs. In the absence of prescriptions for 
determination of a cc and A for particular donor+accretor 
combinations one can consider the "combined" parameter 



a cc A only. We make computations for a cc A=0.5, 1.5, and 
2.5. The choice of a cc A > 1 is also motivated by the fact 
that YLTK considered a formulation of the common enve- 
lope equation which gives for af/ai values that are different 
from the ones given by Eq. J1J and is equivalent to the usage 
of a cc A > 1 in Eq. J1J. Thus, we test both the influence of 
the uncertainties in the definitions of a cc and A and in the 
formulation of the common-envelope equation. 

Generally, we make numerical simulations for the a- 
algorithm, but we also check the 7-algorithm effects on the 
formation of SySs, see T ables ^ and H] For 7, w e consider 
two values - 1.5 and 1.75 (INelemans fc Toutfeoo'jl) . For both 
algorithms we assume that companion mass does not change 
during common envelope phase since this stage is very short. 



2.3 The Model of Symbiotic Stars 

In SySs, the cool component is in FGB or AGB stage. Be- 
cause in the core He burning stage the radius and luminosity 
of the cool component are lower and it can offer very little 
mass by stellar wind for the hot component accretor, this 
stage do es not contribute to the total population of SySs 
(YLTK). We therefore skip this stage in our work. 



2. 3. 1 Mass Loss 

No comprehensive theory of mass loss for AGB stars ex- 
is ts at present. In th is paper, we accept the prescription 
of IHurlev et"aT1 (1200011. In MS Hertzsprung Gap and FGB 
stages, we apply iReimersI il975h mass- loss law 



M(M yr _1 ) = 4 x 10~ 13 



rjLR 

M ' 



(7) 



where r\ = 0.5 and L, R and M are the luminosity, radius and 
mass of the star in solar units, respecti vely. In AGB stage, 
we us e the mass-loss law suggested by Vassiliad is fc Wood! 
(1993) 



log 10 (M) =-11.4+0.0123(P-100max(M/M Q -2.5,0.0)), (8) 

where P is the Mira pulsation period in days given by 

log 10 (P) = -2.07+1.94 log 10 (i?/i? Q ) -0.90 log lo (M/M ). (9) 

When P ^ 500 days, the steady super-wind phase is mod- 
eled by the law 



M(M yr _1 ) = 2.06 x 10" 



,L/L 6 



(10) 



where Voo is the terminal velocity of the super-wind in km 
s _1 ; we use Uoo = 15 km s _1 in this paper. 

If a star loses mass by stellar wind, it loses angular mo- 
mentum too. We assume that the lost matter takes away the 
specific angular momentum of the prospec tive donor. In this 
work , the tidal enhancement of mass loss llTout fc Eggletonl 
1988) is not considered. 

For giant-star accretors we assume that the accretion 
rate is unlimited, since the stars with deep convective en- 
velopes may shrink in the dynamical time scale in response 
to the mass increase. For non-WD accretors accretion rate 
is limited by the thermal time scale: Mi ^ Mi/t^^- For 
larger M2, a common envelope forms. 
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2.3.2 Accretion Rate of Stellar Wind 

Stellar wind accret ion i s crucial for the Sy Ss phenomenon. 
IShima et all lll985|) andlLivio et al l lll98fift showed that for 
SySs the classical iBondi fc Hovld 1)19441) accretion formula 
for the stellar wind is generally valid. The mean accretion 
rate is: 



M hnt = 



-1 ( GMho 



2a 2 (1 + v 2 ) 3 / 



;M CI 



(11) 



where 1 ^ £ w ^ 2 is a parameter (£„ 
is the wind velocity and 



2 

,2 V OTh 2 
V = — V olb 



in this work), v u 



(12) 



where a is the semi-major axis of the ellipse, fJ or b is the 
orbital velocity and total mass M t = Jl4 ot + M coo i. 

If a hot component accretes from the stellar wind of a 
cool component, a fraction of the angular momentum lost 
by the cool componen t is transferred into t he spin of the hot 
component. Following iHurlev et al.l i2002T) . we assume that 
the efficiency of the angular momentum transfer is 1. 



2.3.3 Wind Velocity 

Accretion rate of the stellar wind [Eq. pip], strongly de- 
pends on the wind velocity ti w which is not readily deter- 
mined. Taking into account the fact that the hot component 
may be located in the zone of stellar wind acceleration, we 
assume that the velocity of the stellar wind is related to the 
terminal wind velocity Voo- 

YLTK have defined the wind velocity as 



(13) 



where «oo is the terminal wind velocity and a„ is approxi- 
mated by 



0.04(r/7? d 



1 + 0.04(r/i? d ) 2 



(14) 



where r is th e distance fro m the donor and Rd is the radius 
of the donor. fVoeell fl991) gave an empirical formula for the 
SyS EG Andromedae. This velocity law is: 



Hi 



ci(r/R d ) 



for r/Rd ^ 3.75 



- C2 (r/fl d -c 3 ) for r / Rd > 3 



(15) 



where ci , ci and C3 are 0.2 x (3.75)~ 10 ,§ and 3.42. Below 
we present results of numerical simulations for the values of 
a w given by Eqs. dJ, l|15ft and for a w = 1. 

The definition of terminal velocity in the literature is 
not unique. YLTK take Voo ^CSC which is the surface escape 
velocity. But, as noted bv iHarpeil lll996ft . the main charac- 
teristic of evolved K and early M-giants' cool winds is that 
their terminal velocities are lo wer than the surface escape 
velocity, typically, i>oc ^ Atw. lWinters et"al"l i2000T) found 
two dynamically different regimes for the spherical outflows 
of AGB stars: in the first regime, the terminal wind velocity 
is in excess of 5 km s _1 and the mass-loss rate is higher than 
3 x 10 -7 Moyr -1 ; in the second regime, the terminal wind 
velocity is smaller than 5 k m s -1 and the mass- l oss rate is 
lower than 3 x 10~ 7 M^vr~ 1 . IWinters et all (12003ft fitted the 
relation between mass-loss rates and terminal wind veloci- 
ties derived from their CO (2-1) observations by 



log^MIMoyr" 1 ]) = -7.40 + - log 10 (t 



:>/kn 



(16) 



This rela ti on i s very close to the results of 
lOlofsson et alJ <l2004ft . Mass-loss rates can be obtained 
by Eq . JSJ a nd Vg ^ can be calculated from Eq. \V6\ . 
IWinters et alJ <l2000ft showed the radial structure of hydro- 
dynamic velocity for low and high mass-loss rate models. 
Obviously, Eqs. 1141 or 1151 . which predict the increase of 
outflow velocity, are not adequate for the low mass-loss rate 
model. Semi-regular variable L 2 Pup, whose mass-lo ss rate 
is lower than 3 x 10" 7 Moyr _1 \ Winters et alJl2002ft . is an 
example: its photospheric material is moving at a velocity of 
about 10 km s _1 which is much higher than its terminal wind 
velocity of several km s _1 . In the absence of a clear theory 
of stellar wind generation and behaviour for low-mass giants 
we make numerical simulations for a range of assumptions 
on stellar wind. 

For FGB stars we consider two cases: 

(a) Vac = JjVesc; 

(b) l>oo = V csc . 

For AGB stars we consider cases (a),(b), and case 

(c) in which Voo is determined as follows: 

For mass loss rate higher than 3 x lO _7 M0yr _1 , is deter- 
mined by Eq. 1161 . However, Eq. HIGH is valid for M close to 
10 _6 Mq yr _1 . For higher mass lo ss rate, Eq . 116[l gives too 
high a . Based on the models of lWinters et al\ i2000ft , we 
assume Vac — min(30kms~ 1 , Woo)- Wind velocity is given by 
Eq. 1131 . where a w is defined by Eq. H141 . 
For M ^ 3.0 x lO" 7 M yr~\ we assume that wind velocity 
decelerates from v CBC at the stellar surface to 5 km s _1 at 
r/Rd = 10, using an ad hoc function: 



5km s 1 



^(r/R d ) + ' 



r sC 10.Rd 
r > 107? d - 



(d) In addition, a model with the "standard" terminal wind 
velocity of AGB stars equal to 15 km s is calculated. 



2.3.4 Critical Mass-Accretion Rate of WD 

Most models for eruptions in SySs with WD accre- 
tors involve the nuclear burning of the mat erial ac- 
cre ted by WDs ( Mikola i c wska fc Kenvor] Il992ft. Based 
on iTutukov fc Yungelson! lll976ft and IPaczvriski fc Rudakl 
Jl979ft . SySs powered by hydrogen burning on the surface of 
WDs can be divided further into two subgroups: "ordinary" 
SySs, which are assumed to burn hydrogen steadily if the 
mass-accretion rate is above a critical rate, M st , and sym- 
biotic novae (SyNe) 1 , which experience thermonuclear run- 
aways in their surface hydrogen layers if the mass-accretion 
rate is lower than M st . In SyNe stage and stable hydrogen 
burning stage, the luminosity of the hot component due to 
hydrogen burning is given by the following approximation 
for the core mass-luminosity relation for cold WDs accret- 
ing hydrogen dlben fc TutukovllT996T) : 



L/Lq w 4.6 x lO 4 (M coro /M - 0.26). 



(18) 



For M s t , following YLTK ; we us e the approximation to 
the results of llben fc Tutukovl dl989ft given by 



1 SyNe represents symbiotic nebulae in some literature but sym- 
biotic novae in this paper. 
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log lo M st (M yr- 



-9.31+4.12Mwd- 



1.42(M W d) 2 , 



(19) 



with Mwd in solar units. For steady burning, if the mass- 
accretion rate is above a certain c ritical rate, M cr , an accret- 
ing WD may evolve into a giant ()Keny£n|ll9_81]l) or generate 
an optica lly thick wind f Hachisu ct al. 1996). Based on the 
results of llben fc Tutukovl il989f) . an approximation to M cr 
is given by 

log 10 M cr (M Q yr _1 ) = -9.78 + 9.16M W d - 8.13(M W d) 2 
+2.44(A/wd) 3 , 

where Mwd is in M . If the mass-accretion rate is higher 
than M cr , we assume that the accreted hydrogen burns 
steadily at the surface of the WD and hydrogen-rich matter 
is converted into helium at the rate given by Eq. 12511 be- 
low, while the unprocessed matter is lost from the system 
as an optically thick wind and takes away specific angular 
momentum of the accretor. We also made a test run (case 
13 in Table assuming that excess matter expands and 
transforms the WD into a giant. 

In addition, if the accretion rate is higher than a certain 
value, the release of gravitational energy "mimics" steady 
hydrogen burning. If the luminosity of the WD due to the 
liberated gravitational potential energy, Z/ gra v , is larger than 
101/0, we assumed that the system powered by accretion is 
in the "symbiotic stage" too. I/ gra v (in solar units) is given 
by 



I/grav ~ 3.1 X 10 'Ma, 



Mwd 
Rwu 



ACT 



2 x 10 



-6 ( M v 



(20) 



(21) 



where M acc is in M yr _1 , and Mwd and Rwd are the mass 
and radius of an accreting WD in solar units, respectively. 
We note also that variations in mass-flow rate may cause 
optical outbu rsts that have distinct spectra l features (see 
discussion m iMikoIaiewska fc Kenvonlll992t) . 



2.3.5 Critical Ignition Mass 

For SyNe a certain amount of matter has to be accumu- 
lated prior to the first explosion. Following YLTK we use 
the "constant pressure" expression for AA/ cr it which implies 
that the ignition occurs when the pressure at the base of the 
accreted layer rises to a certain limit: 

-0.8 



(22) 



where Rwo is the radius of zero-temperature degenerate 
objects dNauenberg Ill972l) : 

Rwu = O.O112J? [(M W D/M ch )- 2/3 -(M WD /M ch ) 2/3 ] 1/2 (23) 

with M ch = 1.433M Q and R Q = 7 x 10 10 cm. In fact 
AM^P is a highly sensitive function of the temperature, 
mass, and accretion rate of a WD, as well as of assump- 
tions on the nature of the mixing process at the base of 
the accreted layer and actual nuclear abundances in the ac- 
cretor. Detailed grids of nova mode l s for C O WD-accretors 
may b e found in lPrialnik fc Kovetzl (Il995l) and lYaron et alJ 
(2005). Eq. 1221 agrees with the results of numerical calcu- 
lations for relative ly cold WDs with the temperature 10 7 K 
(lYaron et al l2005lV within a factor of 5, except for the model 
with a 1.4M WD and m = 10 ~ 7 M P) yr" 1 , for which agree- 
ment is within a factor of 7. lNelson et al] i2004h gave numer- 



ical fits to the critical ignition mass es for novae models cal- 
culated bv lPrialnik fc Kovetzl ]l995l) . We made a run of the 
code using the fitting formula for relativel y cold (T ~ 10 7 K) 
WD [see Eq. (Al) of lNelson etafl (120041) in which AM^f 
depends on the mass accretion rates and masses of WD ac- 
cretors]. Its extrapolation to lower mass (0.4 Mq) agrees 
with the results of numerical calculations [by the same code; 
lYaron et al.l JioOS'll to within a factor less than 3. But for 
the lowest mass-accretion rates (10 -12 , lO _12 ' 3 M yr _1 ), its 
errors are larger. Since the systems with such low mass- 
accretion rates provide only a tiny contribution to the total 
number of SySs, we use fits f or 10~ 11 Mfriyr~ 1 in this case. 
Results agree with the ones in lYaron et al] (2005) to within 
a factor of less than 3. 

Because of the complicated dependence of AM cr it on 
the input parameters we ran several simulations varying 
AM crit , see Table [I] 

Oxygen-neon WDs need to accrete a more massive enve- 
lope t4ian_tlie_sarnfijriassCO WDs before the outbursts be- 
gin djose fc Hernanzlll998l) . b ecause of the lower 12 C abun- 
dance in the accretor. Using I.Tose fc Hernanzl (|l99ST ) data, 
we roughly assume that AM cr it for an ONe WD is twice 
that for the same mass CO WD. 



2.3.6 Mass-Accumulation Efficiency 

The efficiency of mass accumulation by a WD can never be 
100%. First, it strongly depends on the strength of symbiotic 
novae and, second, even steady-burning WDs blow stellar 
wind. We define the ratio of the mass of burnt hydrogen 
to the mass of the matter accreted by a WD as oh. The 
strength of a symbiotic nova depends o n the mass and mass - 
accretion rate of the WD. According to lYaron et al] d2005f) . 
we roughly define the boundary between strong SyNe and 
weak SyNe by the mass-accretion rate M ws (in M yr _1 ) 
which is given by 



-11.01 + 6M W d 



logic, M » 



-7.0, 



for Mwd ^ 1M ; (24) 
for Mwd > 1M Q . 



For weak SyNe, i.e., mass-accretion rates being between 
M a t and M ws , we use an a pproximation to qh base d on Fig. 
2 in YLTK [also Fig. 16 of llben fc Tutukovl dl996ty : 



qh 



-4.39- 1.48 log 10 M a , 
-0.10(log 10 M acc ) 2 , 
11.66 + 4.56 log 10 M acc 
+0.45(log 10 M acc ) 2 



for log 10 Macc < -6.36; 



for log 10 Macc > —6.36 



where M aC c is in M yr _1 . 

For strong SyNe, i.e., mass-accretion rates lower than 
M ws , most of the accreted mater is ex pelled and in some 
cases even an erosion of the WD occurs jPrialnik fc Kovet3 
Il995t lYaron et al]l2005^ . Using data on the ejected mass, 
M e j , helium mass fraction in the ejecta, Y e j , and in the con- 
vective envelope, y e nv, heavy element mass fraction in the 
envelope, Z e \ , and in the ejecta, Z e j, given in Table 2 of 
lYaron et al. (2005), the mass of the burnt hydrogen can be 
roughly calculated as: 
for AMcrit ^ M ej , 



AM H = AMcrit xI H - M c j (1 - Y cj - Z ci ) 

-(AMcrit - Mcj) (1 - Kcnv - Z cm ), 



(26) 
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where we assume that the mixing of the accreted layer with 
core material does not occur; 

for AMcrit < M e j, the amount of hydrogen burnt dur- 
ing an outburst is approximated by the difference between 
y^nv and Y e \ which ma y be computed using data from 
iPrialnik fc Kovetd dl995l) . 



AM H = AM CI .it y c 1 



M ci y ci 



(27) 



Then, using data from lYaron et all i2005f) . we fit «h by 



«H = 



AM crit 

= -0.1391 +0.7548Mwd 



(28) 



L{t) = L(Q)t' 



(32) 



where L(0) is given by Eq. (1181 and t is in years. When 
L(t) — IOLq, SySs stage terminates and the cooling time 
tcooi can be obtained. The lifetime of a SyS is the sum of t on 
and t cool . 

We assume that during the outburst WD does not ac- 
crete any matter, due to the high-velocity wind and high 
luminosity. At the onset of decline, WD begins to accrete 
matter from the stellar wind of the cool giant. At the same 
time, the W D loses mass due to intense r adiatively-driven 
stellar wind (lYLTKt liben fc Tutukovll99rj) . We roughly as- 
sume that the rate of the mass loss may be estimated by 



with Mwd in Mq. 

For qh, if the difference between M cr i t and M e j is 
neglected and the greatest part of the energy released in 
the outburst is used to lift the ejected shell from the 
gravitationa l pote ntial well of the WD, then, according to 
lYaron et all 120051. we can give a rough analytical estimate: 



M 



(33) 



qh « 4.54 x 10 



t AfwD 

-Rwd ' 



(29) 



where Mwd is the mass of the WD accretor and -Rwd is 
its radius in solar units, qh in Eq. 1281 agrees with that in 
Eq. J2U within a factor of 3 for O.4M sC M W d < 1.4M 
and the former is usually larger than the latter for a given 
Mwd. 

After a nova, SySs spend certain time in the "plateau" 
stage with high luminosity. This stage lasts for 



t on = 6.9 x 10 



10 Q hAM^P 
L 



yr, 



(30) 



where L is given by Eq. 1181 . Stably burning hydrogen "or- 
dinary" SySs are always in the "plateau" stage, their qh is 
given by Eq. 1251 . In this work the stage of hydrogen burning 
is denoted as "on" -stage. 

For weak SyNe and stable hydrogen burning, we as- 
sume that the burnt hydrogen is deposited at the surface of 
the WD and that the unprocessed matter is lost from the 
system by the stellar wind that takes away specific angular 
momentum of the WD. For strong SyNe, the mass of the 
WD increases if M ei < AM crit . If M ei > AM crit , the WD is 



eroded. We fit the lYaron et alJ i2005T) data by a formula 



i\/„ 



AAf crit 



-1.0372(log 10 M a . 



0.0371 (log 10 M a , 



(31) 



where M acc is in M©yr 1 and which agrees with numerical 
results to within a factor of 1.4. 



2.3.7 Decline of SyNe 

After the outburst, the system remains observable as a SyS 
for a time span tcool, until the WD cools to the temper- 
ature at which its luminosity becomes lower than IOLq. 
iPrialnikl lll986t) sh owed that the luminosit y declines with 
time as I a t~ 114 . Isomers fc Navlorl l|l99fll ) confirmed that 
the obse rved cool i ng rat e of the WD is consistent with the 
model of lPrialriikl 1^)86) . Accordingly, we assume that after 
"plateau" stage the WD enters decline stage during which 
its luminosity decreases as 



where c is the speed of light, v csc is escape velocity and L 
is the luminosity of the WD given by Eq. (l3"2"l . In Eq. $551 . 

v csc is given by 



/2GMwd\ 1/2 

V RwD J 



(34) 



In Eq. 134^ . Rwo is roughly taken as the radius of a cool 
WD. The mass-loss rate given by Eq. 1331 is uncertain 
by a factor of several, since during the main part of the 
"on"-stage the actual radius of the WD is by about an 
order of magnitude larger than the radius of a cool WD 
(see, e.g., iMiirset et al.lll99ll iMurset fc Nussbaumedll994h 
and the performance number of the radiative stellar wind, 
Mv esc /(L/c), is not known. 

Before the luminosity of the WD declines to IOLq, the 
decline stage may be terminated if the amount of accreted 
matter is larger than the critical ignition mass. 



2.3.8 Helium Flashes 

When the helium shell accumulates up to a critical mass, 
the helium flash occurs. This critical mass is about (0.1 
- 0.2) Mq [for initially cold an d non-rotating WDs an d 
M < 3 x lQ- 8 Mp 1 yr- 1 ; e.g., Ilben fc Tutukovl Jl99lfl : 
IWooslev fc Weaverl jl994f) l. After a hydrogen flash, the WD 
becomes "hot" and the helium flashes require ignition mass 
which is much lower than the above-mentioned critical mass. 
For hot white dwarfs we use an approximation for t he criti- 
cal mass sugg ested by Ilben fc Tutukovl Jl989l Il996l) 



AAf«f t = IO^^wSMwd'ACs - 57 , 



(35) 



where M_s is the accretion rate in units of 10 -8 Moyr -1 
and the remaining variables are in solar units. But Eq. 1351 
is not valid for the helium accumulation rate lower than 
lO _8 M0yr _1 . We take the critical mass for helium ignition 
as 0.1 Mq when the helium accumulation rate is lower than 
lO _8 M0yr _1 . Mass accumulation efficiency in the helium 
shell flashes, fo r a wide grid of masses and accretion rates, 
is computed bv lKato fc Hachisul (120041) . We find the growth 
rate of white dwarfs by linear interpolation in these grids. 
We should note, however, that if white-dwarf accretors are 
rapid rotators, the helium accumulation efficiency at low M 
may strongly reduce due t o the dissipation of en ergy at the 
base of the accreted layer jYoon fc Langerll2004l) . 
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2.4 Conversion to Observed Parameters 



n{q) = l, 0<q^l, 



(39) 



To compare our results with observations, we convert lumi- 
nosities of cool components into absolute magnitudes using 
bolometric corrections which can be obtaine d by inter pola- 
tion in the B aSeL-2.0 stellar spectra library of lLeieune et alJ 
Jl997l rT"998l) . SySs are usually detected due to the presence 
of features typical for ionized nebulae in the spectra of cool 
components. This allows us to assume that the simplest se- 
lection effect that governs the observed sample o f SySs i s 
the visual magnitude of the cool co mponents. V r fY TTKf). 
Inspection of the catalogue of SySs feelczyriski et all 2003) 
reveals that the number of objects per unit interval of stellar 
magni tude increases up to V c ~ 1 2.0 mag. In both iKenvorJ 
(1986) and lBelczvhski et al.1 (|2000j), there are about 50 SySs 
with V c ^ 12.0 mag. So, we take 12.0 mag as the limiting 
stellar magnitude of the complete sample and compute the 
number of model objects with V c ^ 12.0 mag. 

The numb er of stars within the distance d from th e Sun 
is taken to be (lEggleton et alJll989t lHan et al.lll995all 



N(d,t) = 



d 3 



(36) 



where t is the evolutionary age of the star, h — 
1200(t/tGai) 1//2 and <GaJ is the Galactic age. The quantity 
do is defined such that 1 / ivdy is the density per unit area in 
the Galactic plane of stars in t he solar neighborhood pro- 
jected onto the plane; following lEggle ton e^ al.| lll98^. we 
take d = 0.054 pc and tn a i=12Gvr. lEggleton et al.l dl989T) 
used the distance d ^ 3 kpc, we roughly extrapolate d to 15 
kpc in this work for comparison with observations. 

To estim ate interst e llar ex tinction, we use the corrected 
relation from Ivan Herkl <196Sft 



A v = 0.14 x cscjfej[l -exp(-10dsin| 



mag.. 



(37) 



where b is the galactic latitude and d the distance in kpc. 
IKenvorJ jl986l) showed that Galactic SySs strongly concen- 
trate toward the plane. We roughly take stellar extinction 
as 1.4 mag per kpc which is the maximum value based on 
Eq. 1571 . 



2.5 Basic Parameters of the Monte Carlo 
Simulation 

For the population synthesis for binary stars, the main input 
model parameters are: (i) the initial mass function (IMF) 
of the primaries; (ii) the mass-ratio distribution of the bi- 
naries; (iii) the distribution of orbital separations; (iv) the 
eccentricity distribution; (v) the metallicity Z of the binary 
systems. 

A simple approximation to the IMF of [Miller fc Scalol 
( 1979) is used. The primary mass is gene rated using the 
formula suggested by lEggleton et al.l 1(1,989) 



Mi 



0.19X 



(1 - X) - 75 + 0.032(1 - X) - 25 ' 



(38) 



where X is a random variable uniformly distributed in the 
range [0,1], and Mi is the primary mass from O.8M0 to 
8M Q . 

The mass-ratio distribution is quite controversial. 
We consider only a constant mass-ratio distribution 
jMazeh et alJll99a iGoldberg fc Mazehlll994l) , 



where q = M 2 /Mi . 

The distribution of separations is given by 

log a = 5X + 1, 



(40) 



where X is a random variable uniformly distributed in the 
range [0,1] and a is in Rq. 

We assume that all binaries have initially circular or- 
bits. We follow the evolution of both components by the 
rapid binary evol ution code, includin g the effect of tides on 
binary evolution dHurlev et alJl2002T) . We take 2 x 10 5 ini- 
tial binary systems for each simulation. Since we present, 
for every simulation, the results of one run of the code, the 
numbers given are subject to Poisson noise. For simulations 
with 2 x 10 5 binaries, the relative errors of the numbers of 
symbiotic systems of different kinds are lower than 7%, with 
exception of SySs with He WD accretors in cases 6, 10 and 
11 (see Table |2l • Under model assumptions for cases 6, 10 
and 11, SySs with He WD accretors hardly form and the 
numbers of them expected in the Galaxy are so small (much 
less than 1 see Tabled that it is impossible to observe them. 
We have made a control run for case 10 with a sample of 
6 x 10 5 binary systems. The maximum difference among the 
results obtained for 2 x 10 5 and 6 x 10 5 samples is less than 
7% of the numbers presented in Tabled except for the re- 
sults on SySs with He WD accretors which agree with the 
numbers given in columns 2, 10, 14 and 22 of Table |21 within 
a factor about 2. However, systems with He-accretors con- 
stitute only a minor fraction of the total SySs population in 
case 10 and, hence, inaccuracy in their number has a neg- 
ligible effect upon our main results. Thus, 2 x 10 5 initial 
binaries appear to be an acceptable sample for our study. 

To calculate the birthrate of SySs, we assume that one 
binary with Mi ^ 0.8Mq is formed annually in the Galaxy 
jYungelson et alJll993l : I Han et al ] |!995aHYLTKh . 



3 RESULTS 

We construct a set of models in which we vary different in- 
put parameters relevant to the symbiotic phenomenon pro- 
duced by hydrogen burning at the surface of WD accretors. 
Table 1 lists all cases considered in the present work, case 
1* is considered as the standard model. In addition to the 
nuclear-burning powered model, we consider an "accretion 
model" that contains systems in which, under the assump- 
tions of case 1, the liberation of gravitational potential en- 
ergy produces symbiotic phenomenon (L grav ^ IOLq) before 
the first outburst of nuclear burning occurs or in the time 
intervals between consecutive nuclear outburst plus decline 
"quasi-cycles" . 



3.1 The history of a typical symbiotic binary 

Our main results concerning the numbers of SySs and oc- 
currence rates of SyNe are given in Table 2. But before we 
discuss them, we present an evolutionary "history" of the 
typical binary that becomes a SyS (Fig. |5J . In this system 
the initial masses of the primary and the secondary are 1.2 
Mq and 1.0 Mq and their initial separation is 1000 Rq. 
We evolve this binary under conditions appropriate to the 
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2500 



Fig. 0i Fig.|2> 

Figure 2. — Luminosity of the primary component vs. evolutionary age in the system with initial masses of components 1.2 and 1.0 
Mq. Panel a represents the evolutionary "history" of the system as a SyS. Panel b shows the details of luminosity evolution during the 
first strong and weak novae experienced by the system (solid and dashed lines, respectively). See text for details. 



Table 1. Parameters of the models of the population of symbiotic 
stars. Case 1* means the standard model. In case 6, Doo of FGB 
stars is taken as v csc and it is ^esc in other cases. In case 7, wind 
velocity ti w is treated as described under item (c) in subsection 
|2 i 3 i 3|_hi_case 13, cri tical ignition mass is taken after Eq. (Al) from 
Nelso n et all <200-fl . 
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on 


case 
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1 


5 


2 ^esc 
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(13) 


AM C WP 


on 
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2 
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on 


case 


4. 
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2 ^esc 
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(13) 




on 


case 


5. 


7 = 1.75 




2 ^esc 


Eq. 


(13) 


A«r 


on 


case 


6. 


a cc \ = 





5 


^csc 


Eq. 


(13) 


AM™ 


on 


case 


7. 


OccA = 





5 


Eq. rm 


Eq. 


(13) 


A«r 


on 


case 


8. 


OiccA = 





5 


15km s 


Eq. 


(13) 




on 


case 


9. 


OiccA = 





5 


2 ^esc 


Eq. 


(14) 




on 


case 


10. 


«c C A = 





5 


2 ^esc 


1 




am wd 


on 


case 


11. 


OiccA = 





5 


2 ^esc 


Eq. 


(13) 


3AM™> 


on 


case 


12. 


a cc A = 





5 


2 ^esc 


Eq. 


(13) 




on 


case 


13. 


« CC A = 





5 


2 ^esc 


Eq. 


(13) 


* am wd 


on 


case 


14. 


OccA = 





5 


2 ^osc 


Eq. 


(13) 


AM WD 


off 



B (in the latter regime no erosion of WD is expected, see 
§ 12.3.61 . The system experiences 3 weak SyNe, before the 
secondary leaves the FGB stage at point C. After return to 
the giant branch but before reaching point D, the secondary 
is an AGB star but the system does not manifest itself as 
SySs. At point D, a strong symbiotic nova occurs. After two 
strong SyNe, the system begins to undergo weak symbiotic 
novae at point E. After 6 weak SyNe, the system enters the 
stage of steady hydrogen-burning by the accretor at point F. 
At point G, the secondary overflows its Roche lobe and the 
system leaves SySs' stage. After the 4th outburst in the D-F 
stage in Fig|2] in the inter-outburst time intervals, the star 
contributes to the accretion model, since its L never drops 
below 10 Lq . 

Figure [^Ja shows in detail the evolution of the luminos- 
ity during the first strong and weak novae experienced by 
the system. At point A (a), the outburst occurs and the lu- 
minosity rapidly increases. From points B (b) to C (c), the 
system is in the "plateau" stage. From points C (c) to D (d) , 
the system is in the decline stage. At point D(d) the lumi- 
nosity of the accretor drops below 10 Lq and we no longer 
consider the system as a SyS, since its accretion luminosity 
(represented by point E(e)) does not satisfy the criterion for 
SySs. Accretion continues until critical ignition mass of hy- 
drogen is accumulated again. Figure [5]d also shows that the 
"on" -stage for the weak novae is much longer than that for 
the strong ones. 



"standard" case 1*. When the system attains point A in 
Fig. UK, the primary has become a CO WD of M = O.54M , 
the secondary is in the FGB stage, the system still is not a 
SyS. At point A, the first strong symbiotic nova occurs and 
the system becomes a SyS. After the first strong symbiotic 
nova, the system leaves the state of being a SyS, while accre- 
tion onto the WD continues until the next symbiotic nova 
occurs. Before reaching point B, the system undergoes 5 
strong SyNe. As the mass accretion rate increases, accretor 
enters the regime of weak thermonuclear runaways at point 



3.2 Birthrate and number of SySs 

First, we discuss the gross properties of the modeled pop- 
ulation of SyS and then proceed to the more detailed com- 
parison of the influence of different assumptions. As Tabled 
shows, the Galactic birthrate of SySs may range from about 
0.035 (case 6) to 0.131 (case 4) yr _1 and it is about 0.076 
yr _1 in the standard model. Their number is from about 
1,200 (case 6) to 15,100 (case 4) and it is close to 4,300 in 
the standard model. The contribution of SySs with He WD 
accretors to the total number of SySs is negligible, due to 
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Table 2. Different models of symbiotic stars population. The first column gives model number according to Tabled Columns 2 to 8 
give the birthrates of SySs with accretors of different kinds (He-, CO- or ONe WD), the rate of SySs formation through evolutionary 
channels I, II, and III (Fig. and total birthrate. The ratio of the number of SySs in cooling stage and their total number is given in 
column 9. Columns 10 to 13 give the occurrence rates of SyNe with accretors of different kinds (He-, CO- or ONe- WD) and total rate, 
in the 13th column the numbers in parentheses mean the rates of weak SyNe. Total Galactic number of SySs with accretors of different 
kinds (He-, CO- or ONe- WD), the number of SySs formed through evolutionary channels I, II, and III and total number are shown 
in columns 14 to 20, respectively. The 21st column gives the total numbers of Galactic SySs with cool components having apparent 
visual magnitude ^ 12.0. Columns 22 to 25 give the numbers for SyNe with accretors of different kinds (He-, CO- or ONe- WD) that 
are currently in the "plateau" stage of an outburst and their total number, in the 25th column the number in parentheses means the 
number of weak SyNe. 



Model 


Birthrate of SySs(yr" i ) 


N CO oling 


Occurrence Rate of SyNe (yr i ) 


He 


CO 


ONe 


I 


II 


III 


total 


Ntotal 


He 


CO 


ONe 


total 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


Standard 


0.001 


0.070 


0.005 


0.019 


0.008 


0.049 


0.076 


0.32 


0.002 


1.4 


2.0 


3.4 (1.5) 


Accretion 


0.0 


0.045 


0.004 


0.001 


0.008 


0.040 


0.049 












case 2 


0.004 


0.080 


0.006 


0.033 


0.008 


0.049 


0.090 


0.38 


0.012 


1.8 


2.3 


4.1 (1.7) 


case 3 


0.006 


0.086 


0.007 


0.041 


0.008 


0.049 


0.098 


0.41 


0.022 


2.1 


2.6 


4.7 (1.8) 


case 4 


0.018 


0.106 


0.006 


0.074 


0.008 


0.049 


0.131 


0.33 


0.160 


6.9 


6.4 


13.5 (6.0) 


case 5 


0.014 


0.099 


0.007 


0.062 


0.008 


0.049 


0.119 


0.34 


0.120 


3.9 


5.1 


9.2 (3.7) 


case 6 


0.0001 


0.032 


0.003 


0.007 


0.005 


0.023 


0.035 


0.46 


0.0001 


1.1 


0.6 


1.7 (0.9) 


case 7 


0.001 


0.121 


0.007 


0.019 


0.008 


0.102 


0.129 


0.62 


0.002 


5.5 


6.5 


11.9 (2.8) 


case 8 


0.001 


0.078 


0.006 


0.014 


0.004 


0.067 


0.085 


0.59 


0.002 


1.8 


4.1 


5.9 (2.4) 


case 9 


0.003 


0.074 


0.005 


0.026 


0.006 


0.050 


0.082 


0.26 


0.007 


1.0 


1.8 


2.8 (1.3) 


case 10 


0.00001 


0.061 


0.005 


0.003 


0.003 


0.060 


0.066 


0.45 


0.00001 


0.9 


2.5 


3.4 (1.9) 


case 11 


0.0002 


0.047 


0.005 


0.010 


0.008 


0.034 


0.052 


0.17 


0.0002 


0.4 


0.8 


1.3 (0.5) 


case 12 


0.004 


0.089 


0.006 


0.027 


0.008 


0.064 


0.099 


0.52 


0.011 


4.3 


4.5 


8.8 (4.0) 


case 13 


0.001 


0.071 


0.005 


0.016 


0.008 


0.053 


0.077 


0.46 


0.001 


4.8 


4.6 


9.5 (5.7) 


case 14 


0.001 


0.070 


0.005 


0.019 


0.008 


0.049 


0.076 


0.37 


0.002 


1.4 


1.8 


3.2 (1.4) 



— continue 
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10 
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4800 


23 


3 
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5 
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20 
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32 


10 
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6 
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1750 


15100 


73 
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1306 


15 


1920 (1858) 


case 5 
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1120 


7180 
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1750 


11200 


59 


304 


1114 


13 


1432(1380) 


case 6 


<1 


1130 


90 


50 
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380 
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4 


< 1 


103 


1 


105 (99) 


case 7 


^1 


8830 


2470 
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1680 


9430 


11300 


74 


< 1 


1242 


9 


1251 (1167) 


case 8 


^1 


2770 


920 


140 


700 


2860 


3700 


23 


< 1 


301 


7 


309 (280) 


case 9 


^4 


3280 


320 


230 


2120 


1260 


3600 


12 


S= 1 


197 


4 


201 (182) 


case 10 


<1 


1600 


360 


30 


340 


1580 


2000 


6 


< 1 


146 


5 


151 (140) 


case 11 


<1 


3220 


240 


100 


1960 


1420 


3500 


13 


< 1 


393 


5 


398 (383) 


case 12 


5 


5360 


680 


600 


2970 


2470 


6000 


32 


S= 1 


411 


4 


416 (395) 


case 13 


=a 


4980 


420 


250 


2960 


2180 


5400 


26 


< 1 


409 


4 


413 (394) 


case 14 




3290 


400 


250 


2030 


1400 


3700 


21 


< 1 


406 


5 


411 (392) 



their large hydrogen ignition mass, except for cases 4 and 
5 in which their contribution is respectively close to 14% 
and 7%; the contribution of SySs with ONe WD accretors 
strongly depends on assumptions and the ranges from about 
7% (case 4) to 25% (case 8); SySs with CO WD accretors 
contribute the main fraction of the population. Observa- 
tional estimates of the total numb er of SySs range f rom sev- 
eral t housands feovarchuklll970ri to about 30,000 llKenvor] 
119941) or even up to 300.000 jMunari fc Renzinilll992l) . de- 
pending on the assumptions on the distance to a typical SyS 

and on observational sel ection. 

In the catalogue of Bel czvhski et alJ i2000Tl . there are 
52 SySs with V c ^ 12.0 mag. In some cases V c is the magni- 
tude in outburst. In our models, the number of SySs with a 
cool component brighter than V c = 12.0 is from 4 (case 6) 
to 74 (case 7). In the standard model, in combined nuclear 



and accretion models, there are in total 56 SySs within this 
magnitude limit, i.e., about 0.5% of all systems! The situa- 
tion is similar for simulations with other assumptions. Thus, 
we may infer that a sample of observed symbiotic stars that 
may be considered as statistically complete, contains only 
fractions of a percent of all Galactic SySs. 

In the nuclear models, from about 38 (case 7) to 83 (case 
11) percent of SySs are in the "on"-state (including stably- 
burning stars and stars in outbursts) and the rest are in the 
cooling stage. When an outburst occurs, the hot component 
spends some time in the "plateau" stage. We assume that 
the time-span of the "plateau" stage is the lifetime of a SyNe. 
In total, the fraction of the SyNe currently in the "on" -stage 
among all SyS (in the outburst and decline stages) varies 
from about 5% (case 9) to 13% (case 7). 

The contribution of SyNe with He WD accretors to the 
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total number of SyNe is about 31% in case 4 and 21% in case 
5 but it is negligible in other cases; due to very short lifetime 
of outbursts in SyNe with ONe WD accretors, their contri- 
bution is also negligible (its range is from approximately 
0.4% in case 7 to 2% in case 8); most SyNe have CO WD 
accretors. In the Galaxy, the model occurrence rate of SyNe 
is between ranges from about 1.3 yr _1 (case 11) to 13.5 yr _1 
(case 4). The contribution of SyNe with ONe WD accretors 
to the total rate is close to 60% in most cases, in extreme 
cases it is 74% (case 10) and 35% (case 6). The contribution 
of SyNe with He WD accretors to the total rate of Novae is 
negligible. Note, the role of systems with ONe accretors is 
rather uncertain because of uncertainty of their progenitors' 
range and may be one of the reasons for the relatively high 
total occurrence rate of SyNe in our models. 



3.3 Parameters 

1. Common Envelope Evolution: a-algorithm — In the 
cases 2 and 3, the parameter a cc A is increased from 0.5 to 1.5 
and 2.5, respectively. The larger a ce A is, the more easily the 
common envelope is ejected, then the orbital period of a bi- 
nary system after common envelope phase is longer, which is 
favourable for the formation of SySs. The occurrence rate of 
SyNe and the total number of SySs in cases 2 and 3 become 
larger than in the standard model. On a ca X depends only 
formation channel I (Fig.0 in which the common envelope is 
encountered. The number of systems produced through this 
channel increases with a ce A since more systems avoid merger 
in common envelopes. The number of SySs and the occur- 
rence rate of SyNe change with increased a ce X by the same 
proportion as the birthrate in channel I, since the outbursts 
are associated predominantly with relatively more massive 
dwarfs which are produced mainly through channels II and 
III. 

2. Common Envelope Evolution: 7-algorithm — In 
case 5, 7 is increased to 1.75 from 1.5 in case 4. Based on 
Eq. (5), the larger 7 is, the smaller the separation of the 
binary is after the common envelope phase as (since more 
orbital angular momentum is taken from the system). The 
birthrate of SySs having undergone channel I in case 5 is 
lower than that in case 4, but 7's effect is weak. 

However, there is great difference between a-algorithm 
cases (1, 2 and 3) and 7-algorithm cases (4 and 5). The 
birthrate of SySs formed through channel I in case 4 is about 
3.9 times of that in the standard model, the birthrate of SySs 
with He WD accretors in case 4 is about 18 times of that 
in the standard model, and the total birthrate of SySs in 
case 4 is 1.7 times higher than that in the standard model. 
The occurrence rate of SyNe and the total number of SySs 
in case 4 are much larger than those in the standard model. 
The reason is that with 7-algorithm post-common-envelope 
systems are wider than with a-algorithm and this facilitates 
symbiotic phenomenon, allowing more stars to evolve fur- 
ther along FGB or AGB before the second RLOF. The 7- 
mechanism is especially favourable for a larger fraction of 
SySs with He- WD accretors among all SyS since it (i) allows 
the progenitors of He-WDs to avoid merger, (ii) hydrogen 
shells accreted by He-WDs are more massive, while their lu- 
minosity is lower, hence, they live as SySs longer. Also, in 
case 4 more massive white dwarfs survive common envelopes 



and this results in a sharp increase in the occurrence rate of 
SyNe. 

3. Parameter Voo — The comparison of cases 1, 6, 
7 and 8 shows that Voo is a key factor for forming SySs. 
For a given mass loss rate, the decrease of Voo increases the 
efficiency of accretion, facilitating SyS formation. This is 
best exemplified by comparison of cases 1 and 6. In case 
7 in early AGB stage, the accretion efficiency of a WD is 
enhanced due to low Voo, but the mass loss rate of a cool 
giant is not high. In late AGB stage, the mass loss rate of 
a cool giant is high, but the accretion efficiency of a WD 
is low due to rapidly rising v^. As a result of these effects, 
the rate of accretion onto WDs stays in the range which is 
favourable for the occurrence of SyNe (see Fig. |H] and Table 
2). 

4. Para meter q w : — In case 9, the wind velocity law 
suggested by IVogell (Il99ll) is adopted. This empirical law 
gives very low velocities in the vicinity of the cool compo- 
nent, which results in more efficient accretion and hence in- 
creases birthrate and number for the closest systems with a 
He- WD, but slightly reduces them for systems with CO- and 
ONe- WDs. However, for a w = 1 in case 10 wind velocity in 
the vicinity of the accretor is high. This strongly decreases 
the birthrates of SySs formed through channels I and II. 
Also, SySs with He- WD accretors are hardly formed. For 
the systems formed through channel III the birthrate and 
number almost do not depend on q w , since for wide systems 
q w is close to 1 for all wind-acceleration laws. 

5. Parameter AM^d — In cases 12 and 11, AMwd, 
the threshold value for "ignition" of symbiotic novae, are, re- 
spectively, I and 3 times of that in the standard model. The 
occurrence rate of SyNe and total number of SySs rapidly 
increase from case 11 to the standard model and to case 
12, but the numbers of SyNe and SySs in the burning stage 
are basically invariable, since the amount of hydrogen avail- 
able for burning is the sam e in all cases. In case 13, Eq. 
(Al) of iNelson et"afl (l2004h is used for AM$g. Compared 
with Eq. I22H . it agrees within a factor close to 2 when the 
mass accretion rate is lO _9 M0yr _1 . However, Eq. (Al) of 
INelson et alJ <|2004T) gives AM^ that are, on the average, 3 
times lower than the ones given by Eq. (I22H when the mass 
accretion rate is lO _8 M0yr _1 and 7 times lower when the 
mass accretion rate is lO _7 M0yr _1 . So the occurrence rate 
of SyNe increases. 

6. Optically thick wind — Comparing cases 1 and 
14, we find that it is practically unimportant for the model 
of the population of SySs whether there is an optically thick 
wind or not. SyNe are not affected, since they occur at low 
M. In case 14, the total number of SySs decreases because 
in the absence of optically thick winds expansion of accre- 
tors and formation of common envelopes in the highest M 
systems becomes unavoidable. The outcome of common en- 
velope is either merger of components or formation of a 
double-degenerate. 

To summarize, we find that the assumed formalism for 
the common envelope evolution, terminal wind velocity and 
AM^D have the strongest effect upon the occurrence rate 
of SyNe and total number of SySs, introducing uncertainty 
up to factors ~ (3 — 4). The wind velocity law affects the oc- 
currence rate of SyNe and the number of SySs within factor 
«S 2. 
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(a) 



(b) Case 4(7=1.5) 

III 
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Cose 7(v„=Eq. (15)) ( d ) 



1.5 3.0 4.5 6.0 7.5 1.5 3.0 4.5 6.0 7.5 

log i^days) log Pj(days) 

Figure 3. — Gray-scale maps of initial primary mass M; vs. ini- 
tial orbital period P; distribution for the progenitors of SySs. The 
gradations of gray-scale correspond to the regions where the num- 
ber density of systems is, respectively, within 1 - 1/2, 1/2 - 1/4, 

1/4 - 1/8, 1/8 - of the maximum of g log p.Q\ og M . 1 an d blank 
regions do not contain any stars. The cases shown in particular 
panels are indicated in their lower right corners. 
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- Case 7{v = Eq.(15)) 
Accretion Model 
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Figure 4. — Number distributions of SySs as a function of orbital 
periods. Panel a shows distributions in total samples for standard 
model, accretion model and cases 4, and 7, while panel b shows 
distribution for a sub-sample with cool components brighter than 
V c =12.0 mag. See t he key to lines in pa nel a. In panel c ob- 
servat ional data from Mikolaiewska (2003) and Bclczvnsk i et all 
120001) is plotted. 



3.4 Properties of SySs 

In this section, we describe potentially observable physi- 
cal quantities of SySs. We compare the standard model 
(a-algorithm for common envelope evolution), case 4 (7- 
algorithm for common envelope evolution), case 7 (low out- 
flow velocity for AGB stars), and the accretion model. 

Figure |3] shows the distributions of the progenitors of 
SySs, in the "initial primary mass - initial orbital period" 
plane. Most SySs descend from (1.0 — 2.5) Mq stars, this 
is an effect of the IMF. The Figure also shows well that 
for a-algorithm for common envelopes, in most SySs, accre- 
tors form from AGB stars. The distribution log Mi — log Pi 
in case 4 (7=1.5) is different from that in other cases, be- 
cause for the 7-algorithm, common envelopes result in rather 
moderate shrinkage of orbits or even in their expansion, al- 
lowing relatively close systems to avoid merger, as already 
noted above. In this case also stars in FGB contribute to the 
formation of SySs, particularly, by formation of He- WD ac- 
cretors and relatively low-mass CO- WD accretors, see also 
Fig. below. The distribution becomes more even in case 7 
because under assumption of low Voo in the systems formed 
through Channel III velocity of the matter passing by WD 
is lower than in the "standard" case and accretion becomes 
more efficient. For the "accretion" model the distribution 
has more narrow range of periods since, in fact, the objects 
of this model, are a subsample of standard model which had 
L > 10 Lq prior to the first outburst or between outbursts. 



3.4.1 Orbital Periods 

Figure 2] shows the distributions of SyS over orbital periods. 
The plot for the total sample, panel a, shows that in case 4 
(7-algorithm) the distribution is wider than in the standard 
case (cv-algorithm) since in the former case there are sys- 
tems which avoided merger and retained or even increased 
their initial separations. In case 7, there are more SySs with 
long orbital periods than in the standard case thanks to the 
systems that form through Channel III and manifest sym- 
biotic phenomenon because of the low wind velocity which 
enhances accretion efficiency. 

The distribution of observed SyS wit h V ^ 12.0 mag, 
over o rbital periods based o n the data from lBelczvriski et alJ 
feOOOl) and iMikoIaiewskal I^OOot) is plotted in Fig. Hfc. It 
may be compared with the distributions for the model sub- 
samples of systems with V ^ 12.0 mag. (Fig.^fc). The over- 
whelming majority of case 1 and case 4 model systems with 
V< 12 in Fig. [IJj are almost uniformly distributed between 
1000 and 6000. However, orbital periods are measured for 
less than 30 systems out of almost 200 known SySs and it 
is hard to measure long orbital periods. The predictions of 
models about orbital period can be hardly verified. However, 
Fig. 2] suggests that SySs with orbital periods shorter than 
200 day are hardly expected, which is in agreement with 
observations. 



12 Lii, Yungelson and Han 



(a) 




(b) 



I I 



r — v 
L _ it 

i . 



(c) 



0.5 



1.5 



Figure 5. — Number distributions of model Galactic SySs as 
a function of the hot component mass. Fig. [Sk is for the total 
samples, Fig.[SJj is for a sub-sample with cool components brighter 
than V c =12.0 mag, and Fig. I5t is plotted based on the data in 
IMikolaiewskal fcOOaft and lBelczvriski et al.l 12000^ . The key to the 
line-styles representing different models is given in the upper right 
corner of Fig. 15k . 




Figure 6. — Gray-scale maps of the model distributions of the 
masses of hot components vs. orbital periods. The gradations of 
gray-scale correspond to the regions where the number density of 
systems is, respectively, within 1 - 1/2, 1/2 - 1/4, 1/4 - 1/8, 1/8 

— of the maximum of — D 9 £? — — , and blank regions do 

not contain a ny stars. 



systems from Mikolaicwsk 
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is based on the data for observed 



3.4-2 Properties of Hot Components 

In Fig. |S] the distributions of the masses of hot components 
are shown. The distribution for case 4 (7=1.5) is different 
from other cases. In case 4, there is a peak between ap- 
proximately 0.35 and 0.5 Mq which represents the distri- 
bution of masses of He WD accretors. Note, the masses of 
these dwarfs are close to the upper limit of the masses of He 
WDs. In other cases, the contribution of SySs with He WD 
accretors is negligible. For all nuclear models, the distribu- 
tions have peaks around 0.6 Mq and a slight enhancement 
by 1.35 M Q . The peak is due to the distribution of CO WD 
masses. The enhancement at the tail of the distribution is 
due to the massive WD accretors (especially ONe WDs) 
that undergo many outbursts because of low hydrogen ig- 
nition mass. The model distribution for the sub-sample for 
SySs with cool components brighter than 12.0 mag. is shown 
in Fig. |5Jd. The distribution of the estimated masses of hot 
components in systems wit h V c ^ 12.0 mag [as given in 
iBelczvnski et alJ fcOOCft and IMikolaiewskal teOOSt ]. is plot- 
ted in Fig. |5k. Though it is still small number statistics, 
we should note that only with the 7-formalism for common 
envelopes is it possible to obtain the prominent fraction of 
low-mass accretors which is present in the observed popula- 
tion. Note the presence of relatively massive WD-accretors 
also in the observed population, though not so prominent 
as in case 7 and not extending to the range of masses of 
ONe WDs. The excess of massive WDs, if real, may be a 
counterargument to the case 4 model. 

Figure |S| presents, in gray scale, the relation between 



the masses of the hot components and orbital periods of 
SyS. The distributions in Figs. [Sh, and c, are split into two 
domains. The scarcely populated left one represents SySs 
formed through channel I; the right domain is populated 
by systems formed through channels II and III. The gap 
between the two domains is absent in case 4 in which sys- 
tems that went through channel I did not experience a sharp 
decrease of orbital separation. Panel c shows that a low ve- 
locity stellar wind favours symbiotic phenomenon in systems 
within a wider range of orbital separations than in the mod- 
els with high velocity. 

The comparison with observed systems shows that our 
models correctly predict that the majority of hot compo- 
nents of SySs must have masses around (0.5-0.6) Mq, though 
the existence of some lower mass (He- WD) accretors cannot 
be excluded. Note that mass estimates typically assume or- 
bital inclinations of i = 90° or a limit to i (see Table 2 in 
IMikolaiewskal teOOStO . hence, the estimates are lower limits. 
Orbital periods of the systems with measured Mhot are also 
consistent with theoretical expectations. At the same time 
our models suggest the existence of a significant number of 
SyS with yet unmeasured periods larger than 1000 day. 

Panel [fjji shows that there is a group of "outliers" - 
systems with relatively short periods and massive accretors. 
Among them is the prominent recurrent nova system T CrB 
(P w 227.6 day, M hot = 1.2±O.2M , M CO oi = 0.7±0.2M Q ). 
Progenitors of such systems must have high initial mass 
primaries and be sufficiently wide to form a massive WD 
and have large initial mass-ratios of components to enable 
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Figure 7. — Gray-scale maps of the distribution of luminosi- 
ties of hot components in SySs vs. their masses. The grada- 
tions of gray-scale correspond to the regions where the num- 
ber density of systems is, respectively, within 1 - 1/2, 1/2 - 



1/4, 1/4 - 1/8, 1/8 -0 of the maximum of 



and 



aiog LSlog M 

blank region s do not contain any stars. Fig. |7p shows M^ot es- 
timates from iMikolaiewskal l2003fl . Solid line shows relation be- 
tw een mass of accreting WD a nd its luminosity 1181 . dashed line 
is iPaczvnski Il970t IUusll97fl) mass-luminosity relation for AGB 
stars L/Lq = 6.0 x 10 4 (M C /M Q - 0.52). 
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Figure 8. — Number distributions of model SySs as a function 
of the mass-accretion rate of the hot components. 



Table 3. An example of an evolutionary scenario leading to for- 
mation of T CrB-like system. Masses are in Mq, P or b in days. 
The values of binary parameters at the onset of respective stage 
are given. 





MS,MS 


AGB, MS 


CE 


WD,MS 


WD,FGB 


CE 


Mi 


7.09 


6.83 


2.19 


1.28 


1.28 


1.28 


M 2 


0.90 


0.90 


1.25 


1.25 


1.25 


1.19 


-f~orb 


7817 


8155 


10622 


277 


277 


263 



sufficient shrinkage of the orbit in common-envelope event. 
With the notation already used above, we present a numer- 
ical example for formation and evolution of a system that 
may belong to this group (in the standard model) in Ta- 
ble The system manifests symbiotic phenomenon in the 
"WD, FGB" -stage. The giant loses only about O.O6M of the 
matter before RLOF and the WD experiences about 1500 
strong Novae (this number is rather uncertain because of the 
uncertainty in AA/^f ). Note, the mass of WD practically 
does not change. 

Figure shows the relations between orbital periods of 
SyS and the luminosities of their hot components. By our 
definition of SySs, there are 3 states of WDs in them, see 
§ 13. H and Fig. [5] Steady-burning systems and SyNe systems 
in the plateau stage obey the mass- luminosity relation 1181 . 
These systems are located in Fig. along the solid line rep- 
resenting Eq. 1181 . Then, when they are in the decline stage, 
they deviate from 1181 downward. Stars in the plateau and 
decline stages form the upper populated area in Fig. |7Ji,b,c. 
Next, as Fig. [5] shows, stars rapidly transit into an "accre- 
tion" stage where they spend time comparable or even longer 
than in the "nuclear-burning" stages. This explains the for- 
mation of the lower populated domain in Fig.|7|and the exis- 
tence of a gap between the two regions. The lower cut-off at 
L/Lq = 10 corresponds to our definition of SySs as systems 
with a larger luminosity WD. The difference in the density of 
systems in different regions of Fig.|7^,b,c may be understood 
as a consequence of different period distributions and accre- 
tion efficiencies in different cases. This picture is reasonably 
consistent with observations plotted in Fig. |7Ji which shows 
that most observed systems are in the decline stage, but, 
probably, observational selection favours detection of the 
brightest systems only 2 . Stars in the long accretion stage are 
not observed because of their low luminosities. Because hot 
components of the observed systems in the mass-luminosity 
plot are located between Iben-Tutukov and Paczyhski-Uus 
curves but not below the latter curve, Fig. |7|I also suggests 
that mass-luminosity relation for hydrogen-accretvng WDs is 
really close to Iben-Tutukov relation 1181 and differs from 
Paczyhski-Uus relation L/Lq = 6 x 10 4 (M C /M Q - 0.52) 
which describes the luminosity of AGB stars with double 
nuclear-burning shells. 

In Fig. |HJ the distributions of mass accretion rates onto 
the hot components M aC c are shown. For all nuclear models 
the distributions have peaks close to 1O _7 M0 yr _1 ; for the 



2 llben fc Tutukovl 1 1996) also suggest that observed symbiotic 
stars are in the stage when their luminosity is provided by residual 
nuclear burning and cooling luminosity of WD. 
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Figure 9. — Number distributions of model SySs as a function of 
the mass of the cool components. Panel a is for the total samples, 
panel b is for a sub-samples with cool compone nts brighter than 
V c =1 2.0 mag and panel c pre sents data from Bclczvnsk i et all 
l l2000l) and lMikolaiewskal 120031) . 
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Figure 10. — Number distributions of model SySs as a function 
of the mass ratio of components. Fig. llOh is for the total samples, 
Fig. 110b is for a sub-sample with cool components brighter than 
Vc =12.0 mag., and Fig,[Tm presents estimates of q in observed 



ag^ 

systems after iBelczvhski et alj i2000l) and Mikolaicwska (2003). 



accretion model the peak is close to 10~ Mq yr~ . The 
first peak and its extension to the larger L correspond to 
the rate of stable hydrogen burning at the surface of typical 
WD accretors, while the second peak corresponds to the 
release of 101/0 by accretion. 

During an outburst, the hot compo nent may sup - 
port an addition a l high -velocity wind [e.g.. iKenvonl |l986); 
llben fc Tutukovl dl996f) j. As noted by the latter authors, 
wind mass-loss may shorten the lifetime of SyS in the "on"- 
state. Another related effect is the formation of a temporary 
"common envelope" by an extended envelope of the white 
dwarf during outburst and associated mass loss. We neglect 
here these two effects. Note, collision of the winds from com- 
ponents is important for the formation of observational fea- 
tures of SySs, but it does not influence their evolution. 



3.4-3 Properties of Cool Components 

The cool component supplies the matter accreted by the hot 
one. The mass of the cool component gradually decreases. 
In Fig. [5^, the distributions of cool components' masses 
are shown. In all cases, the great majority of masses are 
confined to the range (0.8 — 2.0) Mq. The distributions for 
the model sub-samples with cool components brighter than 
Vc =12.0 mag are shown in Fig. I^Jd, the distribution of the 
masse s of observed cool components of SySs with V c ^ 12.0 
mag. (Bel czvhski et alll2000l : iMikolaiew ska 2003) is plotted 
in Fig. |5Jj. Results of our simulations reasonably agree with 
observations. 

Figure HUh shows the distributions of mass ratios of cool 
components to hot components for the total model samples. 
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Figure 11. — Number distributions of model SySs over mass- 
loss rate of the cool component. Fig. Illb is for to tal model sam- 
ples. F ig. lllb presents the observational data from lSeaauist et alJ 
(1993). " 
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Figure 12. — Number distributions of model SySs over the ratios 
of donor radius to the Roche lobe radius of cool component. 



The distributions for the model sub-samples with cool com- 
ponents brighter than V c =12.0 mag. are shown in Fig. llOb . 
The peaks in Fig. llOb for all cases are close to 1.8. Figur/c HUb 
shows the distribution of the obser vational estimates of q fo r 
systems with V c ^ 1 2.0 mag. from iBelczvriski et al.l |2000); 
iMikolaiew ska (2003). The peak in Fig. [TTJ> is located be- 
tween 2 and 4. We do not consider this discrepancy as criti- 
cal. First, it is based on very small number statistics. Second, 
the masses for observed cool components are typically as- 
signed based on their spectral types and evolutionary tracks, 
but for M-giants the masses within the sa me spectral sub- 
type may differ by a factor larger than 2 iDumm fc Schildl 

In Fig. 1111 the distributions of SySs over mass-loss rates 
from the cool components are shown. The peaks in Fig. Illfa 
for the total samples are between 10~ 7 and lO~ 6 Af0yr _1 . 
The observational e stimates of the mass-loss rates after 
ISeaauist et alJ <1993ft are plotted in Fig. II lb . Comparison 
of Figs. II lb and lllb shows that our results agree reasonably 
well with observations. 

Figure 1121 shows the distributions of the cool compo- 
nents over the ratio of their radii to the radii of their Roche 
lobes. In Fig 1121 the distribution for case 7 is different from 
other nuclear cases. In case 7, the maximum is below 0.3, 
but it is between 0.7 and 0.9 in the standard model and case 
4. The main reason for the difference is velocity of the stellar 
wind. In the standard model the rate of mass-loss increases 
with evolutionary lifetime, a w evaluated by Eq. (1141 de- 
creases, thus, the model favours symbiotic stars with donors 
close to the Roche lobe. In case 7, velocity of the stellar wind 
is low during the early AGB stage and this favours symbiotic 
star s with relatively compact donors. 

iMiirset fc Schmidl (1 1999ft found that M-type giants in 
symbiotic binaries obey the relation R ^ h/2 where R is 
the radius of the cool component and Zi is the distance from 
the center of the cool component to the inner Lagrangian 
point L\. Note, it is difficult to compare Fig. 1121 directly 
with the observational data because Sp-M, Sp-R relations 
are not unique for stars in the late stages of their evolution. 
IMiirset fc Schmidl lll999ft applied, for the rad ii of giants, the 
median values frorn ^umm fc Schildl lll998ft . although the 
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Figure 13. Occurrence rate of SyNe (Fig. I13h ^ and weak 

SyNe (Fig. 113b ') as a function of the duration of outbursts. 



real radius may be several times different . However, if the 
pattern found bv IMiirset fc Schmidl (Il999ft is really correct, 
the model of SySs will require a rather low velocity of stellar 
wind in the low mass loss rate phase but a high wind velocity 
in the high mass loss rate phase, whi c h is co nsistent with the 
stellar wind model of lWinters et alJ (l2000ft . 



3.4-4 Symbiotic Novae 

For the nuclear models, the most important property of 
SySs is their thermonuclear runaways. The length of the out- 
bursts depends on the mass and mass-accretion rate of the 
WD-component. As Fig. 113b shows, SyNe last from several 
months to thousands of years. The distributions in Fig. I13h 
can be roughly separated into two domains. In the first do- 
main, below 1 yr, there are peaks around 4 months. The 
first domain corresponds to the strong SyNe. The second 
domain extends from about 1 yr to several hundred (even 
several thousand) yr. As Fig. 113b shows, to this domain con- 
tribute, predominantly, weak SyNe. 

In the standard model for the population of Galactic 
SyNe, the contribution of weak SyNe to the total occurrence 
rate of SyNe is approximately 40%; it ranges from 24% in 
case 7 [uoo=Eq. 1161 1 to 56% in case 10 (a w = 1). If we 
consider the number of systems experiencing thermonuclear 
runaways, it appears that, due to the shorter lifetime of 
strong SyNe, the contribution of systems with weak SyNe 
to the total number of SyNe-systems is larger than 90%. 

Figure ED shows that some outbursts last for more than 
1000 years. These outbursts occur on low-mass He-WDs. 
This is especially prominent in case 4 (7=1.5) in which a 
large number of He WD accretors form. However, these long 
outbursts are more difficult to detect. Obser ved outbursts 
last for months to years or for dozens of years (lKenvorJl 986: 
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lMikolaiewskall2003l) . The outburst of AG Peg lasted for 100 
years. Thus the double-peak structure of the distribution in 
Fig. ll3h is, crudely, consistent with observations. We suggest 
the existence of two populations of SySs - those with short 
outbursts and those with long ones. 

Note that the observed flashes of SyNe with small am- 
plitude (1 to 3 mag) may also be due to the variations of 
mass-transfer rate and/or accretion disk instab ilities and as- 
sociated variations in the nuclear burning rate fbuschil986t 
[ Bisikalo et ai]l2002l : iMitsumoto et alJl2005t ISokoloski et alJ 

F or the occurrence rate of SyNe, llben fc Tutukp-d 
( 1996), taking into account the fact that four outbursts were 
observed within 4 kpc from the Sun in the past 25 years, 
obtained a mean value of v ~ 3 yr" 1 . From the number 
of SyNe's "fossils" , iLewisl ( 1994j) estimated the frequency of 
SyNe as about 0.4 yr -1 , which may be considered as the 
lower limit. Our model estimates are between 1.3 yr -1 and 
13.5 yr -1 . In the standard model, the o ccurrence rate of 
SyNe is close to 3.4 yr -1 , which agrees with llben fc Tutukovl 
( 1996) estimates. Comparing our results to the observational 
estimates, one should note the following: first, not all SyNe 
are observed; second, it is difficult to confirm the event of 
SyNe if it lasts for more than 100 yr. The problem of the in- 
completeness of the sample of Symbiotic Novae was, in fact, 
never explored. 

3.4-5 Symbiotic Stars and Supernovae la 

All models for Supernovae la (SNe la) involve accreting 
WDs. The efficiency of the accumulation of helium shell 
mass is very important if we consider SySs as the poten- 
tial progenitors of SNe la. Helium flashes may occur in our 
models because the critical ignition helium mass is very low 
(See Eq. 135H ) when the mass accumulation rate of helium 
is higher than lO _8 M0yr _1 . The range of their occurrence 
rate is between about 0.0006 yr -1 (case 6) and 0.008yr _1 
(case 4), and the occurrence rate in the standard model is 

clos e to O.OOSyr" 1 

Kenvon et all (Il993l). lYLTKl I Yungelson et. alJ (11996 



and Yungelson et al.l il998l) showed that it is unlikely that 
SySs can produce SN la via accumulation of the Chan- 
drasekhar mass. Our results agree with this conclusion. In 
all cases, the frequency of events in which a CO WD reached 
the Chandrasekhar limit is lower than 10 yr . The rea- 
son is the lack of massive accreting CO WDs and the inef- 
ficient increase of core masses (Fig.EU and even erosion of 
accretors (especially, for massive WD). This is clearly seen 
from Fig. 1141 accumulated helium shell masses do not exceed 
0.1M o . 

The low efficiency of helium accumulation also makes 
symbiotic stars bad candidates for SNe la produced by 
"edge-lit detonation (ELD)" — detonation of core car- 
bon initiated by detonation of accreted helium in WDs 
in the mass range 0.6 — O.9M0, after accr eting 0.15 — 
0.2M fl of He at M ~ lQ-^M^vr" 1 . [e.g.. dLivnd Il990t 
iLimongi fc Torna,mhellT98^ : IWooslev Weaver! Il 994M. ~fn 
all our models, the detonation of accreted helium can not 
occur, that is, there is not a sample in which accumulated 
helium shell mass exceeds O.15M0. We calculated the rate 
of events in which a 0.1 Mq helium shell is accumulated on 
CO WDs with masses larger than O.6M0 in all models. 
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Figure 14. — Occurence rate distributions of SySs as a function 
of the mass changes of CO WD accretors. The bins to the left 
of at x-axis show eroding WD, while to the right of — mass- 
accumulating WD. 



This rate is < 10 yr" . Furthermore, as we mentioned 
above, in the case of fast rotation, outbursts of accreted He 
may occur after accret ion of only several 0.01 Mq of matter 
iYoon fc Langerll2004l) . Then, "Helium Novae" may occur. 



3.5 Comparison with YLTK 

YLTK performed a detailed p opulatio n-synthesis study of 
SySs. In the standard model of lYLTKl the estimate of the 
birthrate of SySs in the Galaxy is 0.076 yr -1 and their 
number is close to 3,300, which is similar to our results in 
the standard model. Certain differences between the models 
stem from the fact that Eq. Q takes into account the sta- 
bilizing effect of a large mass of stellar core upon RLOF for 
stars with convective envelopes that was not considered in 
YLTK. 

The main difference with lYLTrd is the estimate of the 
birthrate of SySs with He WD accretors. In lYLTKl it is about 
0.024 yr" 1 ; for CO WD accretors it is about 0.049yr _1 . Be- 
cause of difference in lifetimes, the numbers of systems with 
different accretors are more similar: 1560 and 1740. In the 
standard model of the present paper, the number of SySs 
with He WD accretors is negligible. There are two main 
reasons for this: 1. F or the same valu e of a ce , the common 
envelope formalism of lWebbinkl lll984h - Eq. (J^J , gives larger 
reduction of component separation than the YLTK model, 
leaving no room for formation of SySs with He WD accre- 
tors. In case 4 (7=1.5), the common envelope model gives 
a smaller reduction of binary separation than Webbink's 
model, then SySs with He WD accretors appear. 2. In both 
models, binaries with He WDs form through common en- 
velopes or RLOF from FGB stars. However, in YLTK the 
upper mass range of the progenitors of He WDs extends to 
2.5 Mq, while in the present study it is taken as 2.0 Mq. 
This difference also allows additional He-accretor systems to 
form in simulations of YLTK. 

Both studies agree that symbiotic stars are not likely 
precursors of SN la. 
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4 CONCLUSION 

We performed a detailed study of the formation of symbi- 
otic systems with WDs as hot components, employing the 
population synthesis approach to the evolution of binaries. 
Several important conclusions can be drawn: 

1. The number of nuclear symbiotic stars in the Galaxy 
may range from about 1,200 to 15,000. These numbers 
are compatible with observational estimates. The model 
birthrate of SySs in the Galaxy is from 0.035 to 0.131 yr -1 . 

2. The estimated occurrence rate of symbiotic Novae 
is between 1.3 yr -1 and 13.5 yr" 1 ; weak SyNe contribute 
to this rate from about 0.5 to 6.0 events per yr. This es- 
timate greatly depends on the critical ignition mass of the 
hydrogen shell AMyyrj arm accretion efficiency. It may be 
compared to the estimate of t he Galactic occurren ce rate of 
classical Novae 41 ± 20 yr -1 (lHatano et alJll997l) . It looks 
likely that the models with high rates of symbiotic novae 
strongly overestimate this rate. This can mean that (i) the 
relative fraction of massive accretors is overestimated be- 
cause of too steep initial-final mass relation; (ii) the values 
of hydrogen-ignition mass AM^J implied by us may be un- 
derestimated; (iii) mass loss by stellar winds from accret- 
ing dwarfs is underestimated or (iv) accretion efficiency is 
overestimated. The latter issue has to be resolved by gas- 
dynamical calculations of mass-flows in symbiotic systems. 
Note also, that tripling of AM^ in case 11 brings occurence 
rate of SyNe to about 1 per yr; having in mind the paucity 
of observed SyNe this low number may be preferred. Igni- 
tion masses given b y Eq. (|22|l are typically hi g her than the 
masse s computed bv lPrialnik fc Kovet j <1995h : lYaron et alJ 
(2005). This may mean that the latter studies also strongly 
underestimate AM^S- 

There may exist two varieties of symbiotic stars — those 
with short outbursts and those with long ones. 

3. The evolution of symbiotic stars is unlikely to lead to 
SNe la via accumulation of Chandrasekhar-mass or to the 
edge-lit detonations. Helium novae may occur in symbiotic 
systems if the energy is efficiently dissipated at the base of 
accumulated He-layer. 

4. The results of the modeling depend on the assumed 
common-envelopes formalism. Within a-formalism, varia- 
tions of combined parameter a cc A do not strongly affect 
model population, since it is dominated by wide systems 
and systems that experience stable RLOF. Increase of a ce \ 
shifts the systems produced via channel I into the ranges of 
parameters of descendants of channels II and III. 

The main effect of applying 7-formalism is the domi- 
nance of systems formed through common envelopes. Com- 
parison with observations, whenever possible (Table 2, Figs. 
4 - 11), suggests that using 7-formalism it is easier to ex- 
plain variety of parameters of observed SySs. This especially 
concerns existence of short-period systems and the tendency 
of SySs for harbouring relatively low-mass accretors. There 
is apparent excess of Novae in 7-formalism models. However, 
this may be also an effect of the underestimate of critical ig- 
nition masses or the overestimate of the number of systems 
with massive accretors. 

5. The wind velocity law affects the occurrence rate of 
SyNe and the number of SySs within a factor of 2. Observa- 
tional data is still not sufficient to provide real restrictions 
on the models. 



6. To summarize, the models for the population of SySs 
mostly depend on the assumed formalism for the common 
envelope evolution and critical mass for hydrogen ignition 
AMJ^. These factors have the strongest effect upon the 
occurrence rate of SyNe and the total number of SySs, in- 
troducing an uncertainty up to a factor of about 4. 
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